1.open pH.inp file in abaqus.
2.open material property, click user define subroutine,
3.input the paramters as, 0.001,  0.1,  3.4, 0.001,  4.3, 0.05 as the properties relating to PROPS(1),PROPS(2),PROPS(3), PROPS(4),PROPS(5),PROPS(6) in the PHsesitive.for file.
4.load the job in abaqus